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We show how the leading physical and Wilson low-energy constants associated with Wilson 
fermions in lattice gauge theory can be determined individually by using spectral information of 
the Wilson Dirac operator with fixed index at finite volume. The methods are demonstrated in 
simulations with leading-order improved Wilson fermions. In addition to the expected suppression 
of the leading term in Wilson chiral perturbation theory we observe a substantial reduction also of 
the higher-order Wilson low-energy constants. 



I. INTRODUCTION 

One of the challenges facing lattice gauge theory sim- 
ulations near the chiral limit is how to disentangle ultra- 
violet lattice artifacts from chiral behavior. As simula- 
tions with Wilson fermions approach the physical point 
of almost massless u and d quarks [H-1H| , this obstruction 
needs to be overcome. The traditional method is to pro- 
ceed through the Symanzik effective continuum field the- 
ory which includes a series expansion in the lattice spac- 
ing a. Associated with each power of a is a set of higher- 
dimensional operators whose transformation properties 
under chiral rotations determine the effective low-energy 
field theory, Wilson chiral perturbation theory. This pro- 
gram has been carried out to high order For a nice 
review, see, e.g., Ref. Q- 

From the numerical point of view, the introduction of 
a whole new set of Wilson low-energy constants in ad- 
dition to the physical low-energy parameters of QCD is, 
however, a challenge. One has to deal with several new 
fitting parameters, and hence the need to fit more observ- 
ables. Moreover, for typical observables in the so-called 
p-regime of Wilson chiral perturbation theory one finds 
that observables tend to depend on particular combina- 
tions of the Wilson low-energy constants, not directly on 
their individual values. 

In this paper we wish to advocate a new strategy to- 
wards determining the leading physical and Wilson low- 
energy constants which is based on the new analytical 
results for spectral correlators of the Wilson Dirac oper- 
ator [H-[l5j. The idea is very simple, and relies on us- 
ing the wealth of detailed information contained in the 
spectral density and spectral correlation functions of the 
Wilson Dirac operator D\y and its Hcrmitian counter- 
part -D5. Here, as in earlier work on how to determine 
QCD low-energy constants based on the Dirac operator 
spectrum [l||, the fact that one has a new finite- volume 



scaling regime plays a crucial role in increasing the ac- 
curacy of the method. In particular, we shall demon- 
strate how the chiral condensate, the quark mass and 
the Wilson low-energy constants Wq, W7 and W% can 
be obtained from the Wilson Dirac spectrum. It is cen- 
tral for this method that one can divide lattice config- 
urations into sectors of fixed index v, as determined by 
the spectral flow [lT}. When combined with more tradi- 
tional approaches from determining combinations of Wil- 
son low-energy constants in the p-regime by measuring, 
e.g., small differences in the charged and neutral pion 
masses [l8|, [l9| or from other combinations of correla- 
tion functions [20] , one has a series of strong consistency 
checks on the obtained values. 

Perturbative results for the effect of 0(a 2 ) terms on 
the spectrum have been compared to numerical simula- 
tions and used to determine Wilson low-energy constants 
in [2l|. There have also been two initial studies of the 
detailed analytical predictions for the Wilson Dirac op- 
erator spectrum [22], [23| , results of which already looked 
promising. Here we explore the effect of introducing a 
simple 0(a) clover improvement [lij], using the coeffi- 
cient given by leading-order weak-coupling perturbation 
theory As we shall demonstrate below, one of the un- 
expected consequences of this simple 0(a) improvement 
is that also 0(a 2 ) low-energy constants of Wilson chiral 
perturbation theory are substantially reduced. Moreover, 
certain asymmetries in the data of Refs. [U [H| that 
clearly could not be explained by Wilson chiral pertur- 
bation theory up to, and including, 0(a 2 ) effects are also 
substantially reduced by this clover improvement. This 
shows that simple clover improvement even reduces the 
terms of order 0(am), where m is the quark mass. This 
is good news for lattice simulations, and a result that 
could not have been anticipated a priori. All numerical 
data in this paper refer to the quenched approximation. 
It would obviously be of interest to have a confirmation 
of our results with the fcrmion determinant included. 
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Our paper is organized as follows. In the next sec- 
tion we briefly review the theoretical setting and intro- 
duce our lattice set-up. In section III we present our new 
numerical results, and explain how different observables 
can be used to probe individual Wilson low-energy coef- 
ficients. We end in section IV with conclusions and an 
outlook for future work. 

II. THE THEORETICAL FRAMEWORK 

The essential input that allows us to compute the ef- 
fects of 0(a 2 ) on the Wilson Dirac operator spectrum 
is the consistency, shown by Sharpe [25| , of the method 
up to contact terms of order am. As in in Refs. @, Q, 
we use the e-regime counting where to ~ e 4 and a ~ e 2 , 
keeping terms up to and including C(e 4 ) (the different 
possible counting rules have been discussed in Ref. [26| ) . 
In the language of finite-size scaling, it means that we 
are considering a regime in which 

mV, zV, a 2 V 

are kept fixed. Here z is a source of the pseudoscalar 
density -075 tp. 

Using the convention of Refs. @, Q, we write the ef- 
fective Lagrangian up to that order as 

C(U) = -(m + z)ZTrU + -(m- z)T,Titf 

- a 2 W 6 [Tr (U + U^)} 2 -a 2 W 7 [Ti (U - U^)] 2 

- a 2 W 8 Tr([/ 2 + C/t 2 ) j (1) 

and we define the partition function in sectors of fixed 
index v by 

Z u N (m,z;a) = f dU Aet" U e~ v c{u] . (2) 

JU(Nf) 

The quenched or partially quenched limits arc defined in 
the usual way. Operationally, the separation into sectors 
of fixed index v as given above corresponds to a sepa- 
ration based on the spectral flow of the Wilson Dirac 
operator 0. This gives a double motivation for follow- 
ing the spectral flow numerically, since it also provides 
us with the real eigenvalues of the Wilson Dirac operator 
Dw (see below). 

We denote the usual Hermitian Wilson Dirac operator 
by As: 

D 5 = 75 (D w + m) . (3) 

The source z introduced above couples to {p^ 5 ip. It is 
useful because the fcrmion determinant of the QCD path 
integral is then, up to a sign, given by 



det( 75 (£V + m) + z) = det(D 5 + z). (4) 

This prompts one to consider the spectral resolvent of 
the Hermitian Wilson Dirac operator D$, 




where are the eigenvalues of D§. From this, the den- 
sity of eigenvalues of D$ follows: 

p 5 (A 5 ) = ilm[G(-A 5 )]|^ 0+ . (6) 

Let us now briefly review some basic facts about the 
eigenvalues of Dw and D§ [9JL many of which follow di- 
rectly from the early paper [271 ] . As is well known, eigen- 
values of Dw are either real or come in complex conjugate 
pairs. The real eigenvalues of Dw play particularly im- 
portant roles, since they provide a definition of gauge field 
topology as the lattice spacing a is taken to zero. Even 
at nonzero lattice spacing these real modes are special. It 
follows from the definition of D$ that its eigenvalues A 5 
are functions of the mass to: A 5 = A 5 (m). Tuning to to a 
value to c at which a zero eigenvalue occurs, A 5 (m c ) = 0, 
is particularly interesting since 

D 5 (j> = D w 4> = - m c (j> . (7) 

This shows that finding a zero of A 5 (to) corresponds to 
identifying a real mode of Dw- It is straightforward to 
see that the argument runs both ways, so that a real 
eigenvalue of Dw also corresponds to a zero of A 5 (to). 
Moreover, in a compact notation, one can show that for 
given eigenstates \j) of real modes of Dw HJUzIj 

dX 5 Am) 

-&r = (il75|j> (8) 

m— m c 

so that the slope at a crossing of the spectral flow is 
given by the chirality of the state As the continuum 
is approached, this chirality goes to ±1. At any finite 
lattice spacing the chirality vanishes identically for all 
nonreal modes [27j ■ 

The chirality of the real modes is also directly related 
to the index v described above since around the physical 
branch 0, 

^sign(j|7 5 |j) = v , (9) 

real 

a relation which displays the topological nature of the 
index: the number of zero-crossings, counted with signs. 
As the lattice spacing a is reduced the number of mul- 
tiple crossings near the physical branch of the spectrum 
goes to zero. On typical configurations in the present 
study, the definition of v given by the spectral flow as 
shown above agreed with the (improved) naive topologi- 
cal charge of the so-called "Boulder method" [28[ , which 
we applied after six HYP smearings to the gauge fields. 
The few cases were the two determinations differed were 
usually associated with multiple crossings, with the last 
one at a fairly large A^ a/ . 

We now remind the reader of a few basic facts about 
the spectrum of the Hermitian Wilson Dirac operator 
[l7l ]. In the continuum, the spectrum of the Dirac op- 
erator D is chiral: for every nonzero eigenvalue iX there 



is a matching eigenvalue — iX. This holds configuration 
by configuration, and is a simple consequence of the 7- 
matrix identity {7 A ',7 5 } = 0. The corresponding Hermi- 
tian Dirac operator 7 5 (_D + m) is also chirally symmetric: 



7 5 (£> + to) 



± \/X 2 + m 2 (j)± 



(10) 



and the pairing is preserved: for every eigenvalue pair 
±i\ of D there is a pair of eigenvalues ±-\/A 2 + to 2 . For 
the zero modes of L>, where there is no chiral pairing, 
there is also no pairing of eigenvalues of 7 5 (-D + ro): these 
eigenvalues become cither +to or —to, depending on their 
chiralities. The zero modes of D are chiral eigenstates. 
For the Wilson Dirac operator Dw these properties are 
violated. In particular, the spectrum of Dw is in the 
complex plane and it is not chirally paired. Likewise, the 
spectrum of D§ = j 5 (Dw + to) is not symmetric config- 
uration by configuration, although the density /9s(A 5 ) is 
symmetric in the v = sector. 

Analytical predictions for the quenched spectrum of 
D5 have been given in Refs. @, Q, although focus there 
was on the simplifying case when only W& is considered 
nonvanishing. General analytical expressions including 
Wq,Wt and W% were also provided in those references. 
The constants Wq , W7 and Wg have fixed signs @, [H, 
Only for W 6 ,W 7 < and W s > is Wilson chiral 
perturbation theory the effective theory of lattice QCD 
with a 75-Hcrmitian Wilson Dirac operator. 

For the clover-improved case studied here, we shall find 
that a^J\Wi\V <C 1. In that case p$ depends only on the 
combination |W 6 | + \W 7 \ rather than on W 6 an d W7 sep- 
arately @. We shall therefore keep Wj = throughout 
the rest of this paper, and show how we can disentangle 
the dependence on We from that on W%. 



III. SIMULATIONS AND NUMERICAL 
RESULTS 

We use two of the gauge field ensembles already consid- 
ered in our previous work 1221 . Both were generated with 
the Iwasaki gauge action [29( and have a lattice size of 
L = 1.5 fm. Some further properties are given in TableUl 
inclu ding the lattice spacing inferred from r^/a values 
from |30l l3l| using r$ = 0.5 fm to set the physical scale, 
and the number of configurations analyzed with \v\ = 0, 
1, and 2. We have applied one HYP smearing |32| to 
the gauge fields before constructing the clover Dirac op- 
erator with tree- level clover coefficient c sw = 1. It is 
known [33[ that the tree-level clover coefficient is fairly 
close to a nonperturbatively improved value when using 
HYP smeared gauge fields. 

We now turn to an analysis of our lattice data. We 
start by numerical results obtained on the 16 4 lattice. 
As a first test of the impact of clover improvement, we 
show in Fig. Q] the density of real modes of Dw in the 
v = 1 sector. On the same plot is shown earlier data 
from Ref. [22| without clover improvement. It is evident 
that clover improvement as expected shifts the position 
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FIG. 1: The effect of clover improvement on the density of the 
real eigenvalues of Dw- The black line is the clover improved 
data and the red line is the data before clover improvement. 
Not only does the peak shift towards the origin, the clover 
term also makes it much more narrow and symmetric. 





size 


a [fm] 


amo 


v = 


\v\ = 1 


\v\ = 2 


2.635 


16 4 


0.093 


-0.03 


1276 


2257 


1518 


2.79 


20 4 


0.075 


-0.027 


1202 


2128 


1408 



TABLE I: Ensemble of pure gauge configurations considered, 
generated with the Iwasaki gauge action. Listed are /3i w , the 
size in lattice units, the lattice spacing, the bare mass in the 
clover Dirac operator and the number of configurations in the 
sectors with M = 0, 1, and 2. 



of the peak towards the origin. Moreover, it is also clear 
that the clover improvement has reduced the width of the 
peak substantially. In the analytic computation within 
Wilson chiral perturbation theory 0, Q this width is de- 
termined by the terms of order a 2 . Hence we conclude 
that the clover term also has a positive effect on the or- 
der a 2 terms. Finally, the fact that the peak has become 
more symmetric must necessarily be due to an improve- 
ment of the terms of order 0(am) ~ e 6 and higher. As 
our analytical results extend only to order e 4 we unfortu- 
nately cannot quantify this suppression of the asymmetry 
further. Overall this seems to show that clover improve- 
ment works better than we had reasons to expect, see 
also [33| . We now turn to the effect of clover improve- 
ment on the spectrum of D§. Since in the small-a limit 
the real eigenvalues of Dw map directly to the peak at 
A 5 = to, we certainly expect improvement here as well. 

The histograms in Fig. [2] give the eigenvalue density of 
the Hcrmitian Wilson Dirac operator in the v = 0, 1 and 
2 sectors as obtained from the 16 4 lattices. Displayed 
also are two fits of the analytic curves obtained from 
Wilson chiral perturbation theory. With the additional 
Wilson low-energy constants there are in total 5 param- 
eters to fit (TV/ a, mTV, W 6 a 2 V, W 7 a 2 V and W 8 a 2 V). 
Fortunately the analytic results of 0, [22| give us a series 
of insights that allow us to address these constants in 
turn. First of all, as mentioned above, the Wilson low- 
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FIG. 2: The spectral density of the Hermitian Wilson Dirac 
operator on the 16 4 lattices. The numerically determined 
eigenvalues have been rescaled with EV/a = 173 for the com- 
parison to the analytic results. The \u\ = 2 data in the lower 
panel allow us to separate the effect of We from that of W&. 

energy constants Wi have fixed signs We < 0, Wj < 
and W 8 > [II. Moreover, [1,0 in the limit 

where ayJWifF <?C 1 there is a factorization of the eigen- 
modes of D§, while those in the index peak have a specific 



dependence on the Wi 's the rest of the eigenva lue density 
is not affected to leading order in ayfW^fF, see Fig. [5] 
below for an illustration. This has the following most 
useful consequences: 

1 ) The eigenvalue density on the opposite side of the 
index peak is almost continuum-like. 

2) The smallest eigenvalues on this side of the origin 
are located very close to — m and thus offer a clean way 
to extract the quark mass parameter. 

3) The value of £ can be estimated by scaling the A 5 
axis and m until the data on this side of the gap match 
the continuum predictions. 

4) The values of the Wi's are finally obtained from 
their effect on the index peak (see below). 

The fits were made, using the above insights, to the 
\v\ = 1 data and have a\J—W§V = 0.1 and W% = (blue 
line) respectively We = and ay/W$V = 0.1 (red line). 
In both cases T,V/a = 173 and mEK = 2.5. 

As can be seen from the middle panel of Fig. [2] the 
\l>\ = 1 data work well both with We alone (blue curve) 
and with W% alone (red curve). 

Having fixed the values of the low-energy constants by 
the \v\ = 1 data we can now check the two correspond- 
ing predictions from Wilson chiral perturbation theory 
against the v = and \v\ = 2 data. In the top panel for 
v = both predictions do well, however, for the \v\ = 2 
data only the prediction from the fit with Wq alone re- 
produces the peak structure of the would-be topological 
modes. 

In Fig.[3]we plot the analytic predictions (the values of 
the low-energy constants still fixed as above) for the real 
eigenvalues of the Wilson Dirac operator against the real 
eigenvalues of the \v\ = 1 and \v\ — 2 sectors. In the top 
panel the red lines display the predictions corresponding 
to the fit with W% alone. The predictions from the fit 
with We alone are given in the lower panel (blue lines). 
Again the \v\ =2 data for the real eigenvalues clearly 
favors the prediction with We alone. 

The reason why We and W% have such a different effect 
on the analytic prediction with \i>\ = 2 is the following: 
The Wg-term in the e-regime of Wilson chiral perturba- 
tion theory corresponds to a Gaussian fluctuating mass 
0. The (5-peak fr om the continuum at m is therefore 
smeared into a Gaussian with an amplitude that simply 
scales with v. In particular, We does not induce eigen- 
value repulsion between the real modes. On the contrary 
it was shown in [?l[l(| that for small but nonzero a 2 WsV 
the density of the real eigenvalues of Dw takes the form 
of the vxv Gaussian unitary ensemble scaled by Aa 2 WsV. 
The effect of W$ is therefore to induce the eigenvalue re- 
pulsion familiar from random matrix theory between the 
real modes. In this way we get a clear distinction between 
the effect of W 6 and W s ■ 

The fact that the real modes of Dyj and the nearly 
topological peak in the spectrum of D 5 both lead to the 
same conclusion is not accidental. In Fig. [4] we plot the 
real modes of Dw for \v\ = 2 together with the distri- 
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FIG. 3: The density of the real eigenvalues of Dw, rescaled 
with EV/a - 173. Clearly the prediction with We, = (blue 
lines lower panel) gives the better description of the \v\ =2 
data. 



bution of the first two positive eigenvalues of £)g appro- 
priately shifted by the bare mass. The essentially per- 
fect match between the two distributions demonstrates 
that the corresponding eigenvectors are almost chiral as 
is expected for |a 2 Wil^| <C 1, see figure 4 of Q and the 
associated discussion. 

Let us now consider the 20 4 data set. In Fig. \E\ the 
eigenvalue density of the Hermitian Wilson Dirac oper- 
ator is displayed for v = (top panel), \v\ = 1 (middle 
panel) and \v\ = 2 (lower panel). Again we have made 
two fits to the \v\ = 1 data one with Ws = (blue line) 
and one with We = (red line). The values in this case 
are ayf—W^V = 0.06 and Ws = (blue line) respec- 
tively We = and ay/WsV = 0.06 (red line) both with 
Ey/a = 220 and mSV = 2.3. 

These values were determined according to the strat- 
egy outlined above. To illustrate the factorization of the 
behavior used we show in Fig. [5] the fit for v = 1 and 
Wg = together with the continuum, a = 0, curve at the 
sa me quar k mass. As is clear the leading-order effect for 
a^/| Wj| V <C 1 is on the index peak only. 



0.05 0.1 0.15 0.2 

a^ and ok -m 

FIG. 4: The density of the real eigenvalues of Dw for \u\ = 2 
on the 16 4 lattice plotted together with the accumulated den- 
sity of the first two positive eigenvalues of D5 on the same 
configurations shifted by the bare mass. The excellent match 
between the two demonstrates that the corresponding eigen- 
vectors are almost chiral. 

The two fits both describe the \ v\ = 1 data well and the 
prediction for \v\ = also works nicely for both fits. For 
\v\ = 2 the analytic predictions again differ in the region 
of the topological modes. Also here it is the prediction 
with We alone that gives the only acceptable fit to the 
= 2 data. This again allows us to conclude that We 
gives the dominant contribution to the discrepancy from 
the continuum. 

Finally, to summarize the effect of clover improvement: 
Before clover improvement [22I ] we found that the val- 
ues of a-y/lWslV where roughly 0.35 for both the 16 4 
and 20 4 lattice s. Here, after clover improvement, we find 
that a -y/ 1 Ws I V is consistent with zero and that ayjWesjy 
takes the values 0.10 and 0.06, respectively. Moreover, 
the asymmetry of the density of the real modes has been 
substantially reduced by the addition of the clover term. 
This indicates an improvement also at order 0(am). 
While this unexpected order 0(a 2 ) and 0(am) improve- 
ment is an analytic challenge to understand, it is clearly 
good news for lattice QCD. 



IV. CONCLUSIONS AND OUTLOOK 

We have shown how the microscopic eigenvalue den- 
sity of the Wilson Dirac operator can be used to deter- 
mine both physical and lattice-artifact low-energy con- 
stants from lattice QCD simulations. We have demon- 
strated this with quenched clover-improved simulations. 
The method relies crucially on our ability to divide the 
lattice configurations into sectors with fixed index. 

While initial numerical tests [22|, [23| had already 
demonstrated the feasibility of computing the micro- 
scopic spectrum of the Wilson Dirac operator and com- 
paring it to analytical predictions, there were aspects of 
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FIG. 5: The spectral density of the Hermitian Wilson Dirac 
operator on the 20 4 lattices. The eigenvalues have been 
rescaled with T,V/a = 220. In blue and red are two fits (see 
the main text) to the \ v\ = 1 data in the middle panel. In the 
top and lower panel the predictions obtained from these fits 
are plotted against the data. 



the measured spectra that displayed a systematic dis- 
agreement with the analytical predictions. In particu- 
lar, the observed asymmetry in the spectrum of the real 



FIG. 6: Comparison of the microscopic eigenvalue density of 
the Hermitian Wilson Dirac operator in the continuum a = 
(black curve) and with a^/WW = iO.06, W s = (blue 
curve). The vertical black line indicates the position of the 
topological 5-peak. To leading order in a ^y\W~\V only this 
peak is affected by the discretization. 



modes of Dw- This asymmetry has been known to be 
present in data from very early on, and yet it is a firm 
prediction of Wilson chiral perturbation theory up to and 
including 0(a 2 ) that the spectrum must be symmetric. 
A quick analysis reveals that terms involving odd pow- 
ers of am in the effective theory are needed to explain 
such an asymmetry. This is beyond present analytical 
predictions. Here we have improved at order a through 
use of the conventional clover term. Since the effective 
theory contains an arbitrary linear shift in the mass term 
anyway, such an improvement is, in this context, not in- 
teresting in itself. But clearly clover improvement has 
an impact on higher-order coefficients as well. We have 
found that clover improvement in fact substantially re- 
duces the higher-order Wilson constants and makes the 
spectrum of real modes much more symmetric, and hence 
in good agreement with analytical predictions. Because 
of the close connection between real modes of Dw and 
the "topological" (threshold) eigenvalues of D5, this im- 
proves the eigenvalue spectrum of -D5 as well. 

Apart from demonstrating substantially better agree- 
ment between theory and numerical data, we have taken 
the opportunity to explain new ways to measure Wil- 
son low-energy constants based on spectral data of the 
Wilson Dirac operator. In particular, we have shown 
how detailed understanding of how the different 0(a 2 ) 
operators in the Wilson chiral Lagrangian influence the 
spectrum can be used to isolate dependencies on the indi- 
vidual low-energy constants. This should be very helpful 
in determining the numerical values of the constants in 
dynamical simulations. Crucial in this context is the sep- 
aration of configurations into sectors of fixed index i>, the 
Wilson analogue of a topological charge. In each sector 
there are definite predictions with which to compare the 
data. 

The present study should be extended to full-scale sim- 
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ulations with dynamical fcrmions. Analytical predictions 
are already available for two light flavors [l2[ . This gives 
entirely new ways to measure the low-energy constants 
of Wilson chiral perturbation theory, that can be com- 
pared to and combined with alternative approaches us- 
ing more conventional space-time dependent observables 
in the p-regime. Although computationally more com- 
plicated, it would also be most interesting to carry out 
analogous studies of the full complex spectrum of D\y- 
The analytical predictions are available in [ll], Ell- In 
particular, the realization of either the Aoki phase [36| 
or the Sharpe-Singleton scenario is closely linked to 
the complex spectrum of Dy/ (l3| . 

Finally, due to the unexpected higher order improve- 



ment caused by the clover term it would be interesting to 
test the effect of clover improvement also in twisted mass 
lattice QCD. As we have demonstrated here, the spectral 
density of the lattice Dirac operator is an efficient tool 
for such an analysis. The necessary analytical predic- 
tions from Wilson chiral perturbation with a maximally 
twisted mass have recently been derived [37T ]. 
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